Field-induced long-range order in the S* = 1 antiferromagnetic chain 
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The quasi-one dimensional 5* = 1 antiferromagnet in magnetic field H is investigated with the exact 
diagonalization of finite chains and the mean field approximation for the interchain interaction. In 
the presence of the single-ion anisotropy D, the full phase diagram in the HT plane is presented for 
H \\ D and H A- D. The shape of the field-induced long-range ordered phase is revealed to be quite 
different between the two cases, as observed in the recent experiment of NDMAP. The estimated 
ratio of the interchain and intrachain couplings of NDMAP {J' / J ~ 10"'') is consistent with the 
neutron scattering measurement. 
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Since Haldane conjectured the presence of the spin 
excitation gap in the one-dimensional integer-5 Heisen- 
berg antiferromagnets, the S* = 1 antiferromagnetic chain 
has attracted a lot of interest. Many theoretical and nu- 
merical works indicated the disordered ground state of 
the system. The mean field approximation for the inter- 
chain interaction combined with the exact diagonaliza- 
tion of finite chains revealed that such a disordered 
ground state should be realized even in the quasi-one- 
dimcnsional S ~ 1 antiferromagnet if the ratio of the 
interchain and intrachain couplings holds zJ' / J ^ 0.05 
where z is the number of adjacent chains. The criterion 
was satisfied for the material Ni(C2H8N2)2N02(C104) 
1^, abbreviated NENP. In fact the Haldane gap had al- 
ready been observed in the temperature dependence of 
the magnetic susceptibility and the high-field magne- 
tization measurements of NENP. |^,|| 

The effective field theory and the size scaling 
analysis with the finite chain calculation indicated 
that the field-induced transition occurs from the disor- 
dered ground state to the gapless Tomonaga-Luttinger 
liquid phase at some critical external field Hd where 
the Haldane gap vanishes. With the interchain inter- 
action, the gapless phase would become the ordered 
phase. Thus the field-induced antiferromagnetic long- 
range order is expected to appear at i?ci in the quasi-one- 
dimensional S — 1 system. NENP, however, didn't ex- 
hibit such a field-induced spontaneous symmetry break- 
ing, because the staggered moment is forced to appear 
even by a smaller external field through the alterna- 
tion of the y-tensors. [p|-pT| Recently the field-induced 
transition to the long-range order was discovered in the 
high-field magnetization measurement of novel quasi-one- 
dimensional S = \ systems Ni(C5Hi4N2)2N3(C104) ||l|] 
and Ni(C5Hi4N2)2N3(PF6) Ij, abbreviated NDMSZ 
and NDMAP, respectively. The magnetic susceptibil- 
ity measurement indicated the presence of the easy-plane 
single-ion anisotropy D for both systems and the value 
was reported as D 0.3 J for NDMAP. The experimental 
phase diagram of both systems in the HT plane exhib- 
ited a quite different shape of the ordered phase between 



the two cases when the external field H is (a)parallel and 
(b)perpendicular to the principal axis of D. This is prob- 
ably because the universality of the symmetry breaking 
is different between the two cases. Thus it would be in- 
teresting to present the theoretical HT phase diagram. 
In this paper, using the exact diagonalization of finite 
chains and the mean field approximation for the inter- 
chain interaction, the phase diagrams for H \\ D and 
H J- D are presented, to compare with the experimental 
ones and give some suggestions for future experiments. 

In order to investigate the magnetization process of 
the quasi-one-dimensional S = 1 antiferromagnet with 
the planar single-ion anisotropy, we consider the Hamil- 
tonian, 

H = J ^ ^ Sj.i ■ Sj.i+i + J' '^'^ ■ 

j i j i 

where J and J' are the intrachain and interchain cou- 
plings, respectively, j is the index of a chain and j') 
means the sum over all the pairs of nearest-neighbor 
chains. We set J — 1 and use the units such as g^iB — 1 
throughout the paper. Only the XF-like anisotropy 

> is investigated, to compare with the experimental 
results of NDMAP. To consider the two cases {a)H || D 
and {h)H ± D, we set a = z and y, respectively. When 
the antiferromagnetic order is possibly induced by the 
external field H, ccy-plane is an easy plane for (a) while 
X-axis is an easy axis for (b) . Thus the universality class 
of the symmetry breaking should be different between 
the two cases; the former corresponds to the XY model 
{U{\)) while the latter the Ising one (Z2). 

In the magnetization process at low temperatures two 
phase transitions are expected to occur at critical fields 
Hci{T) and Hc2{T). At T = the Haldane gap vanishes 
at Hci{0) and the magnetization is saturated at i7c2(0). 
The antiferromagnetic long-range order in the vertical di- 
rection to H should appear between Hci{T) and Hc2{T). 
Based on the mean field approximation for the interchain 
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interaction, the critical field Hd (T) , as well as Hc2 {T) is 
determined by the solution H of the equation 

x:tiiu)iT,H) = ^, (2) 

where Xs^ib) C^' ^^'^ staggered susceptibility of the 
one-dimensional system (J' = 0) and p {=x, y or z) 
specifies the easy component of the antiferromagnetic 
order. To consider the two cases (a) and (b) we set 
— {z,x) and {y,x), respectively. Using the nu- 
merical diagonalization analysis, we obtain the all the 
eigenvalues and eigenvectors of finite chains and calculate 
''^st(iD) (-^' Since the rotational symmetry doesn't re- 
main in any directions for iJ _L £), the available chain 
length is at most L = 8 under the periodic boundary 
condition because of limitation of the memory size. At fi- 
nite temperature, however, the correlation length is much 
smaller than that at T = ~5 to 6 lattice constants 
for iJ = D = 0). Thus even the result for L = 8 is ex- 
pected to indicate some qualitative features of the bulk 
system. In fact the small oscillation in the H dependence 
of Xst(iD)(-^' ^) '^^^ *° finite size effect only appears 
at r < 0.1 for L = 8. Based on obtained HdiT) and 
Hc2{T) from the calculation for L = 8, we discuss on the 
properties of the phase diagram in the i/T-plane in the 
following. 
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FIG. 1. Hci{T) and Hc2{T) plotted versus T for L =8, 6 
and 4. The finite size effect is revealed to underestimate the 
ordered phase. 

To investigate the finite size effect, we show in Fig. |l| 
Hci{T) and Hc2{T) of the isotropic system [D = 0) with 
zJ' = 0.05 determined by the equation (||) for L —8, 6 
and 4 in the TJT-plane. The lower and upper bound- 
aries of the antiferromagnetic ordered phase are given by 
Hci{T) and Hc2{T), respectively. Those lines also give 
the H dependence of the Neel temperature T^. Fig. |l| 
indicates that estimated Tn is higher for larger L, which 
suggests that the finite-chain analysis tends to underes- 
timate the ordered phase. On the other hand, the mean 



field approximation generally overestimates it. Thus on 
the present analysis the finite size effect and the fluctua- 
tion around the mean field are expected to suppress the 
errors due to each other. 
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FIG. 2. Phase diagrams in the HT-plane for {a)H \\ D 
and {h)H _L D with D = 0.3. 

In the following analysis we put D = 0.3, realistic for 
NDMAP. The phase diagrams in the i/T-plane obtained 
from Xst(iD)(-^' ^) ^'^^ L = 8 are shown in Figs. ||(a) 
and (b), for H \\ D and H ± D, respectively. The phase 
boundaries for zJ' =0.05, 0.04 and 0.03 are given in each 
case. Although Hc2{T) exhibits an oscillation due to fi- 
nite size effect, Fig. ||(a) clarifies an interesting feature; 
Hci{T) and Hc2{T) are quite asymmetric for H \\ D. In 
contrast, H dependence of Tn is close to a symmetric bell 
shape for H ± D. According to our calculation for other 
values of D, the difference in the shape of Tn becomes 
clearer for larger D. With increasing Z?, Tn becomes 
more asymmetric ior H \\ D, while closer to a completely 
symmetric bell shape for H J- D. The asymmetric be- 
havior of Tn for II £) is consistent with the H depen- 
dence of the critical exponent rj in the transverse spin 
correlation function {SqS^) ^ (— l)''r~'' of the isotropic 
S = 1 chain at T = 0. In the ground-state magneti- 
zation process of the system the antiferromagnetic spin 
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correlation measured by ry rapidly grows from Hd to the 
maximum while slowly decays toward Hc2 ■ Thus such an 
asymmetric behavior is expected to be characteristic of 
the system with C/(l) symmetry. On the other hand, the 
case (b) with only Z2 symmetry is supposed to be close 
to a classical spin system where the spin wave excitation 
spectrum around m ~ 1 is similar to that around m ~ 0. 

The phase diagrams for || D and H L D {zJ' = 
0.03) are shown together in Fig. |^. They exhibit several 
characteristic behaviors observed in the experiments of 
NDMAP and NDMAZ; (i)At lower temperatures Hci{T) 
is smaller for H \\ D. (ii)At intermediate H Tn is higher 
for H J- D. (iii)The two Tn curves have two intersec- 
tions, although the upper one has not been observed be- 
cause such a strong magnetic field cannot be realized, 
(i) and (ii) imply that the ordered state resulting from 
the U{1) {Z2) symmetry breaking is less (more) sensitive 
to the external magnetic field, while less (more) stable 
against the thermal fluctuation. Our calculation for var- 
ious values of D indicates that the intersection of two Tn 
curves appears only for D > 0. Thus the observed inter- 
section is also an evidence of the positive D for NDMAP 
and NDMAZ. It is consistent with the sign of D deter- 
mined by the susceptibility measurements. [|13|Jl^ 
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FIG. 3. Phase diagrams for 77 || D and _L D with 
zJ' = 0.03 and D = 0.3. 
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FIG. 4. Calculated phase diagrams for zj'^0.010, 0.012 
and 0.014, shown with the experimental results of NDMAP, 
in the cases of (a) 77 || c-axis and (b)77 1. c-axis. Dashed 
curves are guides for the eyes. 

The qualitative consistency of the present results with 
the experiments of NDMAP encourages us to estimate 
the strength of the interchain coupling J' of NDMAP 
by fitting the calculated phase boundaries to the mea- 
sured one. The calculated phase boundaries {Hci{T) and 
Hc2{T)) for zJ' =0.010, 0.012 and 0.014 are shown to- 
gether with the experimental results of NDMAP in Figs. 
§ (a)i? II D and {h)H _L D, where dashed curves are 
guides for the eyes. We set D = 0.3. The temperature is 
normalized as J=30K and the magnetic field is normal- 
ized such that Hci{0) coincides with the experimental 
results. In Fig. ^ (a) it is difficult to fit the calculated 
Tn curve to the experimental one, because the finite size 
effect is still large. Thus we determine zJ' by the agree- 
ment of the maximum value of Tn as zJ' = 0.012 ±0.002. 
On the other hand, the fit of the calculated Tn to the 
measured one is much better for 77 _L Z? in Fig. ^(b). It 
leads to the same conclusion zJ' = 0.012 ± 0.002. How- 
ever, the error of the estimation might be much larger, 
because we neglect the finite size correction and the fluc- 
tuation around the mean held of the interchain couplings. 
Since the number of the adjacent chains is usually more 
than 2, the result zJ'/J - 0.012 leads to J'/J - 10"^ 
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It is consistent with the most dominant interchain cou- 
phng estimated from the recent neutron scattering ex- 
periment. p4[ Using the same mean field approximation 
with the transfer matrix calculation of the chain, a little 
smaller value oi zJ' / J was estimated. It is a reason- 
able difference, because the size correction of the present 
analysis would shift the result to a smaller value. How- 
ever, it is difhcult to determine which result is better, 
because some analyses beyond the mean field approxi- 
mation should modify it into a larger value. Fig. ^ (b) 
suggests that the maximum of Tn of NDMAP could be 
detected with a little larger external field even for if ± _D. 
It would also enable us to perform a more precise esti- 
mation of J' and confirm the difference in the shape of 
Tn between H \\ D imd H L D. 

Recently a sign of the field-induced long-range or- 
der was also observed in the NMR measurement 
[ p6t of another quasi-one-dimensional S' = 1 system 
(CH3)4NNi(N02)3, abbreviated TMNIN. Since J is 
smaller (J ^ 12K), the full magnetization curve has al- 
ready been obtained. Thus the full phase diagram 
in the i?T-plane will possibly be determined in the near 
future. It would clarify more detailed features of the 
field-induced long-range order. 

In summary, based on the mean field approximation for 
the interchain interaction and numerical diagonalization 
of finite chains, we obtained the phase diagram in the 
i/T-plane of the quasi-one-dimensional 5 = 1 antiferro- 
magnet. The results revealed that the shape of the field- 
induced ordered phase is quite different between H \\ D 
and H L D. Fitting the calculated phase boundary to 
the experimental result lead to the estimation of the in- 
terchain couphng of NDMAP as zJ'/J ~ 0.012, which 
well agrees with the estimation by the recent neutron 
scattering measurement {J' /J ^ 10~^). 
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